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Abstract. We investigate the impact of galactic mass loss triggered by ram-pressure stripping of cluster galaxies 
on the evolution of the intra-cluster medium (ICM). We use combined N-body and hydrodynamic simulations 
together with a phenomenological galaxy formation model and a prescription of the effect of ram-pressure stripping 
on the galaxies. We analyze the effect of galaxy - ICM interaction for different model clusters with different masses 
and different merger histories. Our simulations show that ram-pressure stripping can account for ~ 10% of the 
overall observed level of enrichment in the ICM within a radius of 1.3 Mpc. The efficiency of metal ejection of 
cluster galaxies depends at the first few Gyr of the simulation mainly on the cluster mass and is significantly 
increased during major merger events. Additionally we show that ram-pressure stripping is most efficient in the 
center of the galax;y cluster and the level of enrichment drops quite fast at larger radii. We present emission 
weighted metallicity maps of the ICM which can be compared with X-ray observations. The resulting distribution 
of metals in the ICM shows a complex pattern with stripes and plumes of metal rich material. The metallicity 
maps can be used to trace the present and past interactions between the ICM and cluster galaxies. 

Key words. Galaxies; clusters: general - Galaxies: abundances - Galaxies: interactions - Galaxies: ISM - X-ray: 
galaxies: clusters 
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van Gorkom l2Ul)5|l . Hence the metals are a good tracer for 
the present and past interaction processes between galax- 
ies and the ICM. 

The high level of enrichment in the ICM has triggered 
a lot of discussion about the origin of these metals. Several 
groups here investigated the chemical evolution of the 
ICM with extensive numerical simulations. De Lucia et 
al. (20(7^ and Nagashima et al. (2005| investigated this 
problem using N-body simulations combined with semi- 
analytic galaxy formation and evolution. They conclude 
that most of the metals are ejected by supernova driven 
galactic winds mainly from the most massive galaxies and 
they find that there is a mild chemical evolution since z 
= 1. They do not predict the distribution of the metals 
in the ICM. Tornatore et al. ( 2004) use a smooth parti- 
cle hydrodynamic simulation and emphasize the contribu- 
tion of different supernova types to the enrichment but do 
not investigate in detail the influence of different ejection 
mechanisms. Their results indicate that the simulations 
still lack a feedback mechanism which quenches star for- 
mation and transports metals at low redshifts. 

In contrast to these papers our calculations emphasize 
the transport mechanisms of the enriched material and 
we use a shock capturing method for the hydrodynamic 
treatment (see also Schindler et al. 120051 Domainko et al. 
I2()05|l . Furthermore we are able to trace the distribution of 
the enriched material lost by the cluster galaxies. In this 
paper we investigate the impact of the enrichment of the 
ICM by ram-pressure stripping of cluster galaxies. Other 
processes are studied elsewhere (e.g. Kapferer et al. l2005bl 

mm- 

The effect of ram-pressure stripping on individual 
galaxies was studied analytically and numerically by sev- 
eral groups. A first analytical criterion for ram-pressure 
stripping was already given by Gunn & Gott H1972|l . 
According to this criterion galaxies lose material if the 
force due to ram-pressure stripping exceeds the restor- 
ing gravitational force of the galaxy. One result of this 
consideration is the presence of a stripping radius in ram- 
pressure affected galaxies outside which the galactic gravi- 
tational potential is to shallow to prevent the ISM from be- 
ing stripped. Various authors studied this stripping radius 
with respect to acting ram-pressure and galaxy properties 
(e.g. Abadi et al. [TUM Vollmer et al. HUUTl Roediger & 
Henslcr 2005 ) . In general these authors found a reasonable 
agreement between the analytically estimated stripping 
radius and the stripping radius determined by numerical 
simulations. Additionally Roediger & Hensler (WUb) dis- 
cussed the most appropriate analytical description with 
respect to the numerical results. These authors supple- 
mentary adapt the estimate used by Mori & Burkert 
pnnnii which compares the thermal ISM pressure with ram 
pressure. For the case of a constant ram pressure acting 
on the galaxy Roediger & Hensler l|2005 1 found the closest 
coincidence between analytical description and numerical 
simulation for the estimate introduced by Mori & Burkert 
H2000|l . For the case of a time dependent ram pressure like 
it is the case for galaxies moving on non-circular orbits 



through a galaxy cluster, Vollmer et al. H2001|l concluded 
that the Gunn & Gott H1972|l criterion adopted to the 
galactic inclination angle can explain the HI deficiency of 
the galaxies in the Virgo cluster. Further numerical sim- 
ulations for ram-pressure stripping on galaxy scales were 
presented in several papers (e.g. Abadi et al. 119991 Quilis 
et al. l^nUni Mori & Burkert HUM Vollmer et al. I^IMl 
Toniazzo & Schindler l^nUTI and Roediger & Hensler I^UnSjl . 
In contrast to those papers we present simulations showing 
the impact of ram-pressure stripping of all cluster galax- 
ies on the ICM. We obtain simulations on galaxy cluster 
scale and compute the efficiency and time dependence of 
galactic mass loss triggered by ram-pressure stripping. In 
this paper we focus on the effect of ram-pressure strip- 
ping on the metal enrichment of the ICM. The effect on 
the affected galaxies will be shown in van Kampen et al. 
(in prep.). 

This paper is organized as follows: In Section |2] we 
present the numerical method used. In Section |3| we de- 
scribe the general simulation setup and the properties of 
the model clusters. In Section 01 we show our results. And 
in Section [51 we summarize and discuss the results ob- 
tained. 

2. Numerical Method 

We use the same general simulation setup as described in 
Schindler et al. 

2.1. N-body simulations 

Our starting point is a model for the dark matter distri- 
bution within and around the cluster, for which we use 
N-body simulations. The N-body code is that of Barnes 
& Hut l)198(j|l . adapted to include a halo formation recipe 
that prevents over merging on small scales (van Kampen 
et al. 11999!). This paper assumes a ACDM cosmology, with 
57a — 0.7, rtni = 0.3, h = 0.7, and (Tg = 0.93. Initial con- 
ditions are produced using the constrained random field 
method of Hoffman & Ribak H1991|l . as implemented by 
van de Weygaert & Bertschinger H1996(l . Constraints are 
put on the smoothed density field in the form of peak 
heights for a Gaussian window of 6 Mpc, which produce a 
rich cluster in the center of the box for peaks above three 
times the r.ni.s. density fluctuation at that scale. The N- 
body simulations are obtained within a box of 46 Mpc (for 
Ho = 70 km s^^ Mpc^^) with 64^ particles. The softening 
length of the N-body simulation is 20 kpc which is also 
the resolution of the dark matter distribution. The whole 
cluster is described by tens of thousands of dark matter 
particles, which can easily resolve the radial structure of 
the cluster as well as substructures. The halos of the galax- 
ies are represented by groups of particles with a minimum 
of 10 particles with a half-mass radius of the smallest halo 
of 30 kpc. For these galactic halos we assume an isother- 
mal profile. The galactic dark matter halo gravitation is 
treated without any spatial details by simply applying the 
gravitational potential. 
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2.2. Hydrodynamic treatment 

For the hydrodynamic part of the simulation we use a 
shock capturing, grid based scheme (PPM, CoUela & 
Woodward [TUMI Fryxell et al. HM^ . PPM is a higher or- 
der version of the Riemann solver based Godunov method 
(Godunov^"959). This allows us to study all the dynam- 
ical effects in the ICM connected to subcluster mergers 
with high accuracy. To achieve high resolution at the 
cluster center where most of the stripping is expected to 
happen and also to trace infalling gasrich galaxies from 
outside the cluster, we compute our simulation on multi- 
ple refined grids fRuffert .1992^ . Radiative coohng of the 
ICM is treated according to the coohng functions given 
by Sutherland and Dopita fl993"). Metallicity is used as a 
tracer to follow the gas which is lost by the cluster galax- 
ies. Metallicity is handled as a single quantity therefore 
we do not distinguish between different chemical elements 
(like iron and a elements) but assume solar composition. 
Additionally to the bulk motion of the ICM, diffusion 
should also further disperse the stripped gas. In the pre- 
sented paper we are mainly interested in the distribution 
of the enriched material therefore only the effect of dif- 
fusion on heavy nuclei is investigated. Diffusion strongly 
depends on the mass and the charge of the involved nu- 
clei (see Lang ll9iJ^ Chapter 3.1.2 and 3.1.5). The most 
prominent detection of metals in the ICM is connected 
to the presence of iron. Diffusion will spread iron in the 
ICM within the Hubble time on length scales of less than 
100 kpc. Therefore this process is inefficient in further dis- 
tributing iron on galaxy cluster scales. 

2.3. Galaxy formation 

The halo formation recipe that is part of the N-body code 
allows us to easily generate halo merging histories, and 
also the merging histories of the galaxies occupying the 
halos. Simple prescriptions for gas cooling, star formation 
and feedback mechanisms make up our phenomenological 
galaxy formation model, which is an improved version of 
the model of van Kampen et al. H1999(l . Important for the 
context of this paper is the inclusion of a disk model, which 
is similar in nature to that of Mo, Mao & White (,1998) . 
but with the inclusion of an observed threshold for star 
formation in disks (Kcnnicut 119891 Martin & Kennicutt 
1^1) . The adoption of the Mo, Mao & White approxi- 
mations produce a gas and stellar disk with the same ex- 
ponential scale-length, but different cut-off radii and cen- 
tral densities. Our modeling of the metallicity evolution 
is described in two papers: in Kapferer et al. Ij2005b|l the 
coupled circuit of the evolution of the stellar, cold and hot 
gas components is given, whereas the chemical evolution 
of the stellar populations can be found in van Kampen 
et al. H1999|l . The star formation history of each galaxy is 
constrained by its merger history. Halo-halo mergers and 
galaxy-galaxy mergers trigger short bursts of star forma- 
tion. Metals which are produced in disk stars can either be 
added to the cold gas disk or can be ejected directly into 



the hot gas, or into the ICM for cluster galaxies that have 
their hot gas halos stripped. We choose the extreme sce- 
nario to add metals produced in disk stars to the cold gas 
in the disk which can then be stripped off through ram- 
pressure stripping. We do not replenish the stellar disk 
from stellar ejecta along the path of a galaxy: star forma- 
tion is truncated 'forever' if the whole gas disk is stripped, 
and the galaxy ends up as an SO in case the stellar disk 
is prominent enough, and as an elliptical if it is not. In 
reality, some stellar mass loss will occur also after the star 
formation has ended. Intermediate mass stars will eject 
material in form of AGB winds and supernova la ejecta. 

2.4. Ram-pressure stripping 

To calculate the mass loss of cluster galaxies due to ram- 
pressure stripping we follow the classical Gunn & Gott 
(1972) criterion because VoUmer et al. H2U01|I concluded 
that for non circular orbits the Gun & Gott 1)197211 crite- 
rion additionally extended for the inclination angle of the 
galaxy can explain the HI deficiency of galaxies located in 
the Virgo cluster. According to this criterion material is 
lost by a galaxy if the force due to ram-pressure, exceeds 
the restoring gravitational force. The Gunn & Gott (|1972|l 
criterion is given by: 



^'ram > 27rGo-starO'g 

with 



PlCMt^gal 



(1) 



(2) 



^ram IS thc ram-prcssure, G is the gravitational con- 
stant, CTstar is the Stellar surface density, dgas is the sur- 
face density of the gas, picM is the density of the ICM 
and Ugai is the velocity of the galaxy relative to the ICM. 
We calculate the mass loss of the cold ISM disks of spiral 
galaxies. We assume that the galaxies consist of an expo- 
nential stellar disk and an exponential gas disk. Due to 
the phenomenological galaxy formation model (see Sect. 
12.3(1 the scale length of the stellar disk and gas disk are 
equal but the central densities are different: 



M 



ar.gas 



star, gas ^-r/Ro 



(3) 



Mstar is the mass of the stars, Mgas is the mass of 
the gas and Rq is the disk scale length. With these as- 
sumptions we can analytically determine a stripping ra- 
dius i?strip: 



= 0.5 X In 
with 

Rstvip 
i?0 



GMstarMgas 



(4) 



(5) 



Material which is located outside this stripping radius 
is then lost by the galaxy. The validity of this treatment 
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was shown with numerical simulations and with compar- 
isons with Virgo cluster galaxies by Abadi et al. |T999")- 
The amount of cold gas which can be kept by the galaxy 
(M(i?,trip)) is: 



M(i?strip) = A/gas [l-(x+l)e-^] 



(6) 



Consequently the amount of gas Mgtrip, which is 
stripped off the galaxy is then: 



Me 



strip 



(x + 1) e" 



(7) 



Mass loss due to ram-pressure stripping should in prin- 
ciple also be influenced by the inclination angle of the af- 
fected galaxy. The impact of the inclination angle of the 
galaxy on the stripping process was investigated by Quilis 
et al. H2000|l with numerical simulation. Those authors find 
that there is only little dependence of the effectiveness of 
ram-pressure stripping for galaxies which move close to 
face on through the ICM. For example a galaxy moving 
at an inclination of 20° to the direction of motion suffers as 
much stripping as face on encounters (Quilis et al. 2000). 
There are, however, also arguments which suggest a cer- 
tain impact of the inclination angle on the stripping pro- 
cess. The column density of the gas that has to be pushed 
by the momentum transfer from the ICM grows with the 
cosine of the inclination angle. Also the area of ram pres- 
sure which is acting on the unit area of an inclined galaxy 
scales with the same factor. In this work, due to the lack of 
more elaborate considerations and for simplicity, mass loss 
of galaxies moving inclined as opposed to face on through 
the ICM is scaled with the cosine of the inclination angle 
with respect to the mass loss derived for an uninclined 
galaxy. Stripped material is added locally to the ICM and 
is removed from the galaxies. Note that this treatment 
can be regarded as an upper limit of the galactic mass 
loss due to ram-pressure stripping as it was shown that a 
minor fraction of stripped material might be reaccreated 
back by the stripped galaxy (see VoUmer et al. l2001|l . The 
stripping condition is evaluated at every time step of the 
cluster scale hydrodynamic simulation. The mass loss of 
the galaxies is calculated according to the actual proper- 
ties of the galaxies as it is derived by the galaxy evolution 
model (see Sect. l2.3|l . We also assume that the gas disk re- 
mains truncated after a stripping event. This means that 
galaxies are only further stripped when they experience 
increasing ram-pressure. The stripped ISM in the outer 
part of the disks should in principle be replaced by super- 
novae and stellar winds. However observations show that 
star formation is truncated in ram-pressure affected clus- 
ter galaxies (e.g. Koopmann & Kenney 2004). Hence only 
a passively aging stellar population with a significantly re- 
duced mass loss is present in regions of the galactic disk 
where the ISM is stripped off. Although as already men- 
tioned in Sec. 12.31 intermediate mass stars will eject some 
enriched material, we do not take this additional effect 
into account. 



After the stripping radius is reached the mass loss 
continues due to Kelvin - Helmholz instabilities (Mori & 
Burkert'SDIF). Since Roediger & Hensler (2003) argued 
that the contribution to the mass loss caused by Kelvin - 
Helmholz instabilities is small, we do not account for this 
mass loss in our further calculations. 

Observations of nearby spiral galaxies show that abun- 
dance gradients are present in their disks (see e.g. Henry 
& Worthev ll999|l . The slopes of the gradients depend on 
the morphology of the galaxy (Henry & Worthev 
This feature exhibits a certain influence on the enrichment 
of a surrounding medium due to ram-pressure stripping. 
For the case of galactic metallicity gradients which flatten 
in the outer parts of the galaxy we expect only a limited 
effect on the proposed enrichment process. In this case the 
mean abundance of the gas is close to the abundance of 
the bulk of the gas which is located at large galactic radii. 
But this issue certainly demands further analysis. As a 
first approach we do not consider galactic abundance gra- 
dients in our simulations. In future more refined studies 
we will investigate the actual impact of galactic abundance 
gradients on the ICM enrichment. 

2.5. Supersonic galaxies 

In our simulations several galaxies move with a velocity 
relative to the embedding medium which is higher than 
the local sound velocity. In this case of supersonic galaxies 
the classical Gunn & Gott H1972|l approach has to be mod- 
ified. Supersonic galaxies will form bow shocks in front of 
them and the stripping criterion has to be applied to the 
conditions behind these bow shocks. Hydrodynamic quan- 
tities behind the bow shocks can be calculated as functions 
of the quantities in front of the bow shock according to 
the Rankine-Hugoniot conditions for all Mach numbers 
A^ > 1. Using the equation of continuity at the bow shock 
we can calculate the ratio between densities and velocities 
in front of and behind the bow shock: 



P2 



V2_ 

Vl 



(8) 



Here p are the densities, v are the velocities and sub- 
script I stands for quantities in front of the shock and 
subscript 2 stands for quantities behind the shock. For 
the ram-pressure in front of and behind the shock this 
leads to: 



D 2 

Jiam,l,2 — Pi, 2^1 2 



(9) 



As the velocities behind the shock are smaller than 
the velocities in front of the shock also the ram-pressure 
will be smaller behind the shock due to its dependence 
on the square of the velocity. This means that for super- 
sonic galaxies ram-pressure stripping is less efficient than 
the classical Gunn & Gott P1972') approach would suggest. 
We calculate the Rankine-Hugoniot condition for the ac- 
tual Mach number of each individual galaxy. Therefore 
this treatment is also correct for supersonic galaxy mo- 
tion in galaxy clusters with low Mach numbers. The effect 
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is most pronounced for very fast galaxies. In the case of 
Wgai 3> Cs with Vgai being the velocity of the galaxy relative 
to the ICM and Cg the local sound velocity and treating 
the ICM as an ideal gas with 7 = 5/3 with 7 being the 
ratio of specific heats, the ram pressure is by a factor of 
4 lower behind the galactic bow shock than in front of 
the bow shock. An additional effect in supersonic ram- 
pressure stripping could come from the fact that the bow 
shock is bent according to the Mach number. Since our 
simulations are performed on galaxy cluster scale we can- 
not resolve individual galaxies and the form of their bow 
shocks. But it is expected that for low Mach numbers the 
opening angle of the bow shock is quite large. Therefore 
we do not take into account the effect of the bending of the 
bow shock. It is also important to mention that Roediger 
& Hensler l)2(J(J5|) found a better agreement between ana- 
lytic estimates and numerical simulations of ram-pressure 
stripping if a Rankine-Hugoniot correction is applied. 



dependent on the environment it is expected that the to- 
tal amount of material stripped from the affected clus- 
ter galaxies will depend on the properties of the cluster. 
Galaxy clusters with higher masses strip galaxies more ef- 
ficiently than galaxy clusters with lower masses. This is 
due to a higher ICM density and to a higher velocity dis- 
persion in high mass systems. 

In our simulation we derive the total amount of metals 
stripped from the cluster galaxies. By comparing the total 
mass of metals stripped we test the dependence of galactic 
metal loss on cluster mass. Indeed we find that the total 
mass loss is highest in the high mass cluster (model cluster 
1) and lower in the low mass clusters. We see that in the 
high mass cluster in our sample (model cluster 1) the total 
amount of metals ejected into the ICM is about four times 
higher than in the model clusters with lower mass. 



3. Simulation set up and model clusters 

We run the hydrodynamic simulation on four levels of 
nested grids centered on the cluster centers, the largest 
grid being (20 Mpc)^ and the smallest (2.5 Mpc)^. Each 
of the nested grids has a resolution of 128^ grid cells. We 
start the hydrodynamic simulation at a rcdshift of z = 1 
for numerical reasons and follow the evolution of the ICM 
to a redshift of z = 0. This is the time interval of interest 
as ram-pressure stripping is expected to be most efficient 
after the galaxy clusters have formed. 

In order to test the impact of ram-pressure stripping on 
different galaxy clusters we choose model clusters which 
span a wide range of masses and evolutionary histories. 
Our model clusters contain no cD galaxies and do not 
feature cooling cores at their centers. The model clusters 
are as follows: 

Model cluster 1: This is the most massive cluster of our 
sample. The final total cluster mass at z = inside a 
radius of 3 Mpc is 1.3 x lO^^M©. This cluster experi- 
ences only small merger events after z = 1. 

Model cluster 2: This is the least massive cluster in our 
sample. The final total cluster mass at z = inside a 
radius of 3 Mpc is 7.4x lO^^M©. This cluster undergoes 
a major merger with a mass ratio of 1:3 at a redshift of 
0.5. Additionally it features the infall of several smaller 
subsystems. 

Model cluster 3: This cluster has a mass of 8.7 x 
IO^'^Mq inside a radius of 3 Mpc at the final redshift 
of z = 0. It undergoes two minor mergers. 

4. Results 

4.1. Total amount of stripped metals 

The first quantity of interest is the total amount of ma- 
terial and metals which is ejected into the ICM by ram- 
pressure stripping. As ram-pressure stripping is a process 



0.002 - 




[Gyr] 



Fig. 1. Temporal metallicity evolution Afmctais/A/icM in 
the different model clusters. The values are given for all 
galaxies which are located inside the central (2.5 Mpc)'^ of 
our simulations. For an explanation of the different slopes 
for the different model clusters see Sect. 14.2! 

Another interesting quantity is the fraction of the ICM 
which actually originates from cluster galaxies and is lost 
by them due to ram-pressure stripping. To investigate this 
question, we compare the mass of the hot ICM, which in 
our simulations has never been inside any galaxy, to the 
amount of material which has actually been lost by the 
galaxies (see Figl^J. 

We find that about 10% of the ICM in the central (2.5 
Mpc)^ of our simulation originates from ram-pressure af- 
fected galaxies. The fraction of the ICM which originates 
from cluster galaxies is comparable for all our model clus- 
ters (see Fig. l4.2|) . This result is expected since in the high 
mass cluster also the mass of the ICM is larger than in 
the clusters with lower mass. Mixing then the ICM in the 
high mass cluster with a larger amount of material which 
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Fig. 2. Fraction of the ICM which originates from cluster 
galaxies and is lost by them due to ram-pressure stripping. 
The values are given for the central (2.5 Mpc)^ of our 
simulations. For an explanation of the different slopes for 
the different model clusters see Sect. 14.21 

is lost by the galaxies results in a comparable fraction of 
lost ICM to primordial ICM as in the low mass clusters. 

4.2. Mass loss rate due to ram-pressure stripping 

We now investigate the combined mass loss rate of all 
cluster galaxies as a function of time. It is interesting to 
mention that the time evolution of the mass and metal 
ejection is different for different model clusters (see Fig. 
n] and Fig. ^ . With the analysis of the total mass loss 
rate of all cluster galaxies as a function of time we further 
explore this behavior. The mass loss of all cluster galaxies 
as a function of time is shown in Fig.|2|and Fig2] 

We find that in the first 3 Gyr the mass loss rate in 
the high mass cluster (model cluster 1) is about an order 
of magnitude higher than the mass loss rate in low mass 
clusters. There is also a prominent peak at the start of the 
simulation in model cluster 1. This first peak originates 
from the mass loss of gas rich galaxies which experience 
the effect of ram-pressure for the first time by switching on 
the calculation. After 3 Gyr the mass loss of the galaxies 
in the high mass cluster drops significantly. We interpret 
this as an effect of cluster mass. Stripping in the high mass 
cluster can reduce the size of the gas disks of galaxies quite 
fast and then this mechanism becomes less efficient (see 
Fig. El). 

The situation is very different for the merging clus- 
ter (model cluster 2). In the case of a cluster experienc- 
ing mergers the prehistory of ram-pressure is also followed 
within the infalling subsystems. In this system mass loss 
starts at a quite small rate but between 3 Gyr and 5 Gyr 
the rate rises by an order of magnitude. In this time in- 
terval also the major cluster merger happens. Here we see 
the influence of the merger history on the evolution of the 



Fig. 3. Combined mass loss rate of all galaxies as a func- 
tion of time. The values are given for the central (2.5 
Mpc)^ of our simulations. In this plot the influence of 
cluster mass on the combined stripping rate of all clus- 
ter galaxies can be seen. Model cluster 1 is about twice as 
massive as model cluster 3. 



Model cluster 2 

Model cluster 3 




time [Gyr] 



Fig. 4. Combined mass loss rate of all galaxies as a func- 
tion of time. The values are given for the central (2.5 
Mpc)^ of our simulations. In this plot the influence of the 
merger histories of the galaxy clusters on the combined 
stripping rate can be seen. Model cluster 2 undergoes a 
major merger at t — A Gyr. 

mass loss rate due to ram-pressure stripping (see Fig^. At 
about 6 Gyr several galaxy groups fall towards the cluster 
center and the mass loss rate increases again. Subcluster 
merger can raise the stripping rate by an enhanced ve- 
locity dispersion, increased ICM density and due to infall 
of gas rich galaxies onto the main cluster. The strongly 
enhanced mass loss rate during subcluster mergers is also 
the reason for the quite different slopes for model cluster 
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1 and 3 and model cluster 2 in Fig. ^ and Fig. |21 The 
influence of cluster mergers on ram-pressure stripping of 
cluster galaxies is further analyzed in Mair et al. (in prep.). 

Finally we note that model cluster 3 features a simi- 
lar mass loss rate of its galaxies (with some scatter due 
to small merger events) over the entire simulation time of 
8 Gyr. Hence on cluster scale ram-pressure stripping is a 
process which can act over a quite long period of time. 
This supports the idea of some enrichment over the whole 
redshift range of cluster evolution by ram-pressure strip- 
ping (Finoguenov et al. OP 00,) . 

4.3. Contribution to the chemical enrichment of the 
ICM 

Gas originating from cluster galaxies is enriched with 
heavy elements due to feedback processes from stars and 
supernovae. If a large amount of this gas is lost by the 
galaxies, the ICM into which this gas is mixed will also be 
enriched up to a certain level. In this section we investigate 
the contribution of ram-pressure stripping on the chemical 
evolution of the ICM. Tornatore et al.(f2DDT) have argued 
that some galactic mass loss and metal transport at low 
redshifts is needed to explain the overall chemical abun- 
dance of the ICM. Ram-pressure stripping is a promis- 
ing candidate which could account for the necessary mass 
transfer away from the cluster galaxies. 

We explore the capability of ram-pressure stripping for 
enrichment at low redshifts. The mean observed chemical 
abundance of the ICM is about 0.3 solar. For compari- 
son with the observed values we derive the mean chemical 
abundance in the central (2.5 Mpc)^ of our simulation. 
The final level of enrichment is 0.03 solar for model clus- 
ter 1, 0.05 solar for model cluster 2 and 0.04 solar for 
model cluster 3. This corresponds to a contribution of 10 
- 15% to the observed overall enrichment of the ICM, de- 
pending on cluster mass and evolutionary history of the 
cluster. In case of the galaxy cluster which experiences a 
major merger (Model cluster 2), the enrichment happens 
later during the cluster evolution in comparison to the 
non-merging cluster (e. g. Model cluster 1). Here it is in- 
teresting to note that the fraction Mmctais/MjcM (see Fig. 
nj is rising faster then the fraction Mgtrip/MicM (see Fig. 
12) for model cluster 2. This effect results from the inter- 
nal metallicity evolution of the cluster galaxies. Material 
which is stripped later from galaxies is enriched to a higher 
level since stellar feedback is acting for a longer period of 
time on this gas. 

4.4. Metallicity profiles 

For a better comparison of the results from our simula- 
tions with X-ray observations we derive emission weighted 
metallicity profiles. These are shown for the three model 
clusters in Fig. |5| The profiles for all three model clusters 
are remarkably similar and do not show a strong signature 
of the different merger histories of the different systems. 
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Fig. 5. Emission weighted metallicity profiles of the three 
model clusters at a redshift of z = 0. Ram-pressure strip- 
ping can account for the observed metallicity in the central 
100 kpc in all simulations, but at larger radii other pro- 
cesses than ram-pressure stripping have to contribute to 
the enrichment of the ICM. 

In the calculated profiles we clearly see the highest 
level of chemical enrichment of the ICM at the cluster 
center. At the cluster center the environmental impact on 
the galaxies will be most important so it is expected that 
there the enrichment is most pronounced. In the central 
100 kpc ram-pressure stripping can account for an enrich- 
ment of about 0.2 solar which is close to the observed 
value. At larger radii the chemical abundance drops very 
quickly and other processes than ram-pressure stripping 
must contribute to the chemical evolution of the ICM. 

4.5. Distribution of enriched material 

Metallicity can be used as a tracer to identify the location 
of the gas which was lost by cluster galaxies. To investigate 
the distribution of the enriched material and to compare 
our simulations with observed galaxy clusters, we obtain 
emission weighted metallicity maps. For a discussion on 
the effect of diffusion on the metallicity distribution see 
Sect. 12.21 In contrast to the full 3D distribution of the 
metal rich areas, enriched material located in high den- 
sity regions is much more pronounced in emission weighted 
metallicity maps. The emission weighted metallicity maps 
clearly show a complex pattern of the chemically enriched 
material (see Fig. A nonuniform distribution of the 
enriched material is in good agreement with recent obser- 
vations (e.g. Hayakawa et al. I2()()4|l . 

In all the model clusters we see stripes of metal rich 
material which were left behind by infalling galaxies. Some 
of them are even several hundred kpc long. It is interesting 
to note the similarity of these stripes with the plume of 
stripped gas observed at the center of the Virgo cluster 
in HI (Oosterloo & van Gorkom l2005|l and with the trails 
of ionized gas behind two galaxies in Abell 1367 (Gavazzi 
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Cluster 1 t = 4.0 Gyr Cluster 1 t = 7.9 Gyr 




Cluster 3 t = d.O Gyr Cluster 3 t = 7.9 Gyr 




Fig. 6. Emission weighted metallicity maps of the model cluster. The size of the boxes is 5 Mpc on a side. The level 
of enrichment is given in solar units. 
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et al. I2UI)1|I . In contrast to these observations, material 
in our computed enriched stripes is already mixed with 
the hot ICM. Evaporation of stripped gas clouds at the 
center of our model cluster (Tjcm ~ 10® K, njcM ~ lO"*^ 
cm~^) can occur within 10^ yr (VoUmer et al. l2(JUl|l . which 
is the typical time step of the hydrodynamic simulation. 
Stripped gas from cluster galaxies is indeed seen in the 
hot (T = 10® K) ICM of the Coma cluster (Finoguenov et 
al. I2()()4|l in X-ray observations. 

The non-merging model cluster 1 shows metal maps 
with an almost spherically symmetric appearance. There 
is only little evolution between 4.0 Gyr and 7.9 Gyr in our 
simulation (see Fig. two panels at the top) . This is in 
good agreement with the results of Sect ion In contrast 
model cluster 2 shows a very complex distribution of the 
enriched material and there is a lot of enrichment hap- 
pening between a simulation time of 4.0 Gyr and 7.9 Gyr. 
The subcluster merger can clearly be seen in the panel 
showing the situation at 4.0 Gyr. The emission weighted 
metallicity maps for model cluster 2 are presented in Fig. 
121 in the two central panels. Several stripes of enriched ma- 
terial can be seen in the maps of model cluster 3 (see Fig. 
El lower two panels). There is also some evolution up to 
the end of the simulation at 7.9 Gyr. 

Summarizing, our simulations enable us to derive 
abundance maps of the ICM. With the help of these maps 
we can identify the origin of highly enriched regions which 
have recently been observed with X-ray observations (e.g. 
Hayakawa et al. 12004(1 . Furthermore we can investigate the 
location and the fate of the stripped material and compare 
this to up to date X-ray data (e.g. Finoguenov et al. i2004i) . 

5. Summary and discussion 

In the present paper we investigate the chemical enrich- 
ment of the ICM due to ram-pressure stripping of cluster 
galaxies from z = 1 to z = 0. We find that ram-pressure 
stripping plays a significant role in this stage of cluster 
evolution. In particular we investigate the efficiency, time 
dependence and spatial distribution of this process: 

— More than 10% of the ICM originates from gas of clus- 
ter galaxies which was lost by them by ram-pressure 
stripping. 

— Enrichment by ram-pressure stripping contributes to 
the overall metal enrichment of the ICM by about 10 
- 15% within a radius of 1.3 Mpc. It is more difficult 
to enrich high mass clusters with high ICM densities 
than low mass clusters with lower ICM densities. This 
effect is hardly compensated by the higher efhciency 
of ram-pressure stripping in high mass systems. 

— The total mass loss rate of all cluster galaxies depends 
on the cluster mass. High mass clusters can reduce the 
gas disks of their galaxies very quickly and enrichment 
of the ICM by ram-pressure stripping can take place 
within a few Gyr. 

— Enrichment by ram-pressure stripping persists over the 
entire simulation time and can therefore account for 



some late enrichment of the ICM. This effect is more 
pronounced in low mass systems. 

— Cluster mergers have an important impact on enrich- 
ment by ram-pressure stripping. Subcluster mergers 
can raise the total mass loss of all cluster galaxies sig- 
nificantly 

— Ram-pressure stripping can account for the total ob- 
served level of enrichment within the central 100 kpc 
but at larger radii other processes have to contribute 
to the chemical evolution of the ICM. 

— Metallicity gradients are not very well suited to inves- 
tigate the complex distribution of metals in the ICM. 
Metal maps are much more useful to investigate the 
present and past interactions between the ICM and 
the cluster galaxies. 

— We find an inhomogeneous distribution of the enriched 
material in the ICM. This is in agreement with recent 
observations (e.g. Hayakawa et al. I2004|l 

— In our simulations we see stripes of metal rich material 
which were left behind by ram-pressure affected galax- 
ies. We note the remarkable similarity of the shape 
of these stripes with the plume of stripped material 
observed in the Virgo cluster (see Oosterloo & van 
GorkomUnnSI). 

It seems that merger activity of the galaxy cluster has 
an important impact on the enrichment of the ICM due 
to ram-pressure stripping of cluster galaxies. Studying a 
larger sample of merging clusters will lead to a better 
understanding of this problem. In the future we plan to 
investigate the effect of various merger scenarios on the 
chemical evolution of the ICM. 

In forthcoming more refined studies we will also inves- 
tigate the impact of additional effects like abundance gra- 
dients in the galaxies, stripping due to Kelvin - Helmholz 
instabilities, shapes of bow shocks of supersonic galaxies 
and different galaxy evolution scenarios on the enrichment 
of the ICM due to ram-pressure stripping. 
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